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Abstract 



The retarded Van der Waals force between a polarizable particle and a per- 
fectly conducting plate is re-examined. The expression for this force given 
by Casimir and Polder represents a mean force, but there are large fluctua- 
tions around this mean value on short time scales which are of the same order 
of magnitude as the mean force itself. However, these fluctuations occur on 
time scales which are typically of the order of the light travel time between 
the atom and the plate. As a consequence, they will not be observed in an 
experiment which measures the force averaged over a much longer time. In 
the large time limit, the magnitude of the mean squared velocity of a test 
particle due to this fiuctuating Van der Waals force approaches a constant, 
and is similar to a Brownian motion of a test particle in an thermal bath 
with an effective temperature. However the fluctuations are not isotropic in 
this case, and the shift in the mean square velocity components can even be 
negative. We intèrpret this negative shift to correspond to a reduction in the 
velocity spread of a wavepacket. The force fluctuations discussed in this pa- 
per are special case of the more general problem of stress tensor fluctuations. 
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These are of interest in a variety of areas fo physics, including gravity theory. 
Thus the effects of Van der Waals force fluctuations serve as a useful model 
for better understanding quantum effects in gravity theory. 

PACS numbers:34.20.-b, 05.40. +j, 42.50Lc, 12.20.-m 
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I. INTRODUCTION 



The retarded Van der Waals forces between pairs of àtoms and between an atom and 
a perfectly conducting plate were first calculated by Casimir and Polder [1]. In the long 
distance limit, where the àtoms may be described by a static polarizability, these forces may 
be interpreted as due to shifts in the vacuum energy of the quantized electromagnètic field. 
This is most clearly illustrated by the Casimir effect [[|, which may be viewed either as the 
retarded Van der Waals force between a pair of perfectly conducting plates, or as the shift 
in the vacuum energy due to the plates. It has recently been measured accurately || ||. 
Similarly, the Casimir-Polder force between a plate and an atom has been confirmed by a 
experiment by Sukenik et al 0. Note that the large distance limit of the theory can be 
applied to any polarizable particle, but not just an atom. 

Because these forces have their origin in the vacuum fluctuations of the electromagnètic 
field, it is perhaps not surprising that the forces themselves are fluctuating forces. The first 
discussion of the force fluctuations was given by Barton f?]||, who considered fluctuations 
of the Casimir force between plates. In this approach, one considers a spatial and/or time 
average of the force. It is found that the fluctuations diverge in the limit that the averaging 
time goes to zero. Further work along the same lines was done by Eberlein ||. Jaekel 
and Reynaud |ÏC| have also discussed Casimir force fluctuations, especially for accelerating 
mirrors, using an approach based upon fluctuation-dissipation theorems. In this paper, 
we will consider the fluctuations of the force between an atom and a perfectly conducting 
plate from an approach somewhat different from that adopted by either of the above sets of 
authors. Our approach is based upon the Langevin equation. The solution of this equation 
to find the mean squared velocity of the particle involves a time integration which introduces 
a natural averaging scale. We will show that this averaging is suflicient to yield finite results. 

The problem addressed in the present paper can be viewed as a special case of the 
larger problem of understanding the quantum fluctuations of the stress tensor [p^ -15], This 



problem is of interest for a variety of reasons, ranging from radiation pressure noise in an 
interferometerJTH], to quantum fluctuations of spacetime geometry driven by stress tensor 
fluctuations |T3,16|. 

This paper is organized as follows: The Van der Waals force is reviewed in Sec. |I| and 
then the force -force correlation function will be calculated in Sec. |ITT[ In Sec. [IV], we use 
this correlation function to study the velocity fluctuations of a test particle. Here it will be 
useful to use a decomposition of the correlation function into three parts, and to study the 
effect of each part individually. Our results will be summarized and discussed in Sect. [V]. 



II. THE VAN DER WAALS FORCE 

First, let us recali the result for the mean force. We assume that the atom can be 
described as a point particle with a static polarizability a. Its interaction energy with a 
classical electromagnètic field, E, is 

í/ = --E 2 . (1) 

We will use Lorentz-Heaviside units with h — c — 1, but will restore factors of h and c in key 
results. We now assume that the electromagnètic field is quantized, and that its quantum 
state is such that (E) = 0. However, (E 2 ) ^ 0, and there is a mean force given by 

(F) = |V(E 2 ). (2) 
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Quantities such as (E 2 ) in the presence of a plate may be calculated from the photon 
Hadamard function: 



1 



= G^(x, x') = ~ (A^x) A v {x') + A v {x') A,(x)) = G$ + G 



' IXU' 1 



where 



Att 2 D(x,x') 

is the Hadamard function for empty space, and 

4ir 2 D(x } x') 



Gfxvl 



(3) 



(4) 



(5) 



is an "image" term due to the presence of the conducting boundary jll|,|ï2| . Here rj^ = 
diag(—l, 1, 1, 1) is the Minkowski mètric, and z is the unit vector in the z direction 

z" = (0,0, 0,1). (6) 

Furthermore, D(x,x') is the squared geodèsic distance between x and x', 

D{Xj X ') = _ (í _ ff + ( X _ X 'f + ( y _ y'f + ( Z _ Z 'f j (7 ) 

and -D(x, x') is the corresponding distance between a; and the image point x! = (f, x', y', —z'), 
D{x, x') =-{t- t'f + {x- x'f + {y- y'f + (z + z'f . (8) 

The vacuum expectation value of a product of elèctric fields in the presence of a con- 
ducting plate is given by 



{Ei Eji) — do <9q' Giji + di dy Goo 1 — do dy Gío> — di do> Goj>- 



(9) 



We are using a notation in which unprimed indices refer to the spacetime point x, and 
primed indices to the point x'. Thus, do = d/dt, dy = d/dx^ , etc. The quantity (EiEy) 
is divergent in the coincidence limit, x' — > x, but the divergent part does not contribute to 
the force in Eq. (^|). For the calculation of the mean force, the only part which is of interest 

is the renormalized expectation value, which is obtained when G^i rather than is used 
in Eq. This is simply subtracting out the pure vacuum contribution, and is same as 
normal ordering with respect to the Minkowski vacuum. Equation (|9|) can be rewritten as 



(Ei Ef) = (: Ei Eji :) + (Ei Ey)o , 
where the normal-ordered term is 

(: Ei Ej/ :) = do do' G^y + di dy Goo> — do dy G^ — di do> Goj' 
and the vacuum term is 



(10) 
(11) 



(Ei Ey)o = d d > Gfy + di dj> Gfçj, — d dy Gf QI — di d ' G °y . (12) 

If we combine Eqs. @, and (|ÏT1) , we obtain the Casimir-Polder result for the mean 
force: 
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This is an attractive force in the direction perpendicular to the conducting plate. 

Recali that Eq. ( f[3|) strictly holds only when the particle is described by a static 
(frequency-independent) polarizability. For the case of a one electron atom in its ground 
state, Casimir and Polder gave a more complicated expression which reduces to Eq. (|T3| ) in 
the large z limit. In the case of a macroscopic particle with nontrivial dispersive properties, 
there is the possibility of having a force which is either attractive or repulsive, and larger in 
magnitude than that given by the above expression ||17|| . In the present paper, we will deal 
only with the case of a frequency-independent polarizability. 



III. THE FORCE-FORCE CORRELATION FUNCTION 

Now we wish to study the fluctuations in this force. This may be done by examining the 
correlation function (: F(x) :: F(x') :) and the expectation value of the squared force (: F : 2 ). 
However, we will encounter the quantity (: E 2 (x) :: E 2 (a/) :), which is formally divergent 
in the coincident limit x' — > x. Unlike the quadratic expectation vàlues encountered in 
the case of the mean force, we cannot simply render this quantity finite by subtracting 
its expectation value in the Minkowski vacuum state. Following the method used in our 
previous works [HJT5I, this two point function can be decomposed into three different terms 



by using Wick's theorem 

(: E 2 (x) :: E 2 (x') :) = (: E 2 (x)E 2 (x') :) + (: E 2 (x) :: E 2 (x') :) cross + (: E 2 (x) :: E 2 (x') :)„ , 

(14) 

which are the fully normal-ordered term, the cross term and the pure vacuum term, respec- 
tively. In the coincidence limit x — *> x f , the fully normal-ordered term is a well-defmed local 
quantity. The cross term contains a state-dependent divergence, but can be made finite 
with careful regularization in the integral. The pure vacuum term is also divergent in the 
coincidence limit, but it is state- independent and cancels when we measure the difference 
due to a changes of the boundary condition. The fully normal-ordered term can be expressed 
explicitly as 

(: E 2 (x)E 2 (x') :) = (: E % E l E y E 3 ' :) 

= (: Ei E l :) (: E f E j ' :) + (: E t E r :) (: E i E j ' :) + (: E t :) (: E? E r :) (15) 

and the cross term is 

(: E 2 (x) :: E 2 (x') :) cross = (: E0 : : E f E j ' :} cross = 4(: E t E f :)(^") . (16) 



The physical content of both of these terms has been discussed by us [|ï4] , p!5[| in other contexts. 
In general, both terms can contribute to the fluctuations of the stress tensor, or other 
quadratic operators. 

The force- force correlation function (: Ei :: Fy :), evaluated at x = x' but at different 
times, can be obtained by the formula 



a 2 



(: F k {t, z) :: F k (t', z) :) = — d k 8 k > (: E^ :: E f E> :) (17) 
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Again, this correlation function contains two parts we are interested in, namely the fully 
normal-ordered term and the cross term, and the cross term is divergent in the coincident 
limit, t — > t'. The contribution from these two terms will be examined separately in the 
following section. The idea is to investigate the velocity dispersion of a test particle due to 
this fluctuating Van der Waals force. 



IV. VELOCITY FLUCTUATIONS OF A TEST PARTICLE 

We can better understand the effects of these fluctuations by studying the motion of 
partides subjected to the fluctuating force, which will be described by a Langevin equation. 
Consider partides which start at rest at time t = 0. The mean velocity at a later time t is 
given by 

(:v k (t) :) = i jf* : > > ( 18 ) 

where k — x,y, z. To simplfy the analysis, we assume that the distance of the particle from 
the plate does not change significantly in a time t, so that z is approximately constant. Then 
the dispersion around the mean velocity in the fc-direction at a later time is given by 



Aüfc(í)) = \ f dh f dt 2 \(: F k (t u z) : : F k (t 2 , z) :) - (: F k (t u z) :) (: F k (t 2 , z) :}1 . (19) 
m z Jo Jo 1 J 



rt rt 

dti \ 
Jo Jo 

This can be decomposed into two terms 

(A«2> = (: A^:) + (A^) cross , (20) 
which are the fully normal-ordered term 

(: Av 2 h (t) :) = \ [ dt x [ dt 2 [(: F k {t u z) F k (t 2 , z) :) - (: F k (t u z) :) (: F k (t 2 , z) :)1 (21) 
m z Jo Jo 1 J 

and the cross term 

(A^(t)) cross = f dt x I dt 2 (: F k (h, z) : : F k (t 2 , z) :) cross . (22) 
m z jo jo 

Here the pure vacuum term is dropped because we are only interested in the difference due 
to a change of boundary conditions, which is the change caused by adding a plate. The fully 
normal-ordered term and the cross term will now be discussed in turn. 



1. The fully normal-ordered term 



Consider the force fluctuations due to the fully normal-ordered term, Eq. (|I5[ ). The 
diagonal components of the force-force correlation functions are 



[:F k (t,z)F k (t',z) :) 



2 r 



d k d ¥ (: E i E i EjiE j ' :) 



(23) 



Note that the off-diagonal terms will be zero in the limit, x — > x'. Use Eq. (j|) and Eq. (|ÏT|). 
and we find that the elèctric field two point functions can be expressed as 
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(■ E X E X > :) 
(: E y Eyi :) 



2tt 2 
1 

2^ 



2 

D 2 

2 

D 2 



4(x - x'f 4(í - í') 2 



D 3 
I) 3 



+ 



D 3 
4(í - t') 2 

D 3 



and 



•1 



2vr 2 



2 

D 2 



Plug Eqs. (ED 

(: F x (í, z) F a (t', ^) :) = (: F„(í, z) F y (í', z) :) 
and 



4(z + z') 2 + 4(t-í') 21 



D 3 D 3 
D and © into Eqs. (0) and (0). We find 

4a 2 (5T 4 + 16T 2 z 2 + 48 z 4 ) 



(24) 
(25) 

(26) 



7T 4 (T 2 — 4z 2 ) ' 



(27) 



(: F z (t,z)F z (t',z) :) - (: F z (t,z) :)<: F 2 (t',z) :) 



4a 2 (5 T 6 + 308 T 4 z 2 + 944 T 2 z 4 + 1728 , 



vr 4 (T 2 - 4z 2 Y 



where T = t — t' and the product of the mean force is 



:)(:F,(í',z) :) = - 



(28) 



(29) 



All of these results are independent of x and y, and are Lorentz invariant under boosts in 
the directions parallel to the plate. In the limit t' — > t, these fluctuations become 



(: AF x 2 (z,í) :) = (: AF*(z,t) :) 



3« 2 



256 vr 4 z 10 



and 



27 a 2 



256 vr 4 z 



4 r 10 



(30) 



(31) 



Here the z component is about 5 times the x and y components. If we compare the expecta- 
tion value of the squared force from Eq. (31) with the square of the expectation value from 
Eq. (|ï"3|), we obtain a measure of the force fluctuations: 



(: F z (t } z) 2 :) - (: F z (t,z) :> 2 



(:F z (t,z) :Y 



3 
4 



(32) 



which is of order of unity and shows that the force is fluctuating considerably. Note that 
even though there are no mean forces in x and y directions, the deviation of the force in these 
directions are still non-zero. Furthermore, the correlation function, Eq. ( p7[ ) and Eq. (p8"D, 
becomes small if T ^> z, i.e. for time separations large compared to the distance of the atom 
from the plate. This shows that the characteristic fluctuation time is of the order of z. 

However, the behavior at time scales larger than the characteristic fluctuation time z is 
also important, and is needed to find the velocity fluctuations. Note that the force correlation 
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functions, Eqs. (|27|) and (0), are singular at T = 2z, a time separation equal to the round- 
trip time light travel between the particle and the plate. This singularity is presumably an 
artifact of our assumption of a perfectly reflecting plate, and would hence be smeared out in 
a more realistic treatment. However, we will see that the integrals can be made well-defmed 
even with this singular integrand. 

Let us now change integration variables to T = t\ — t 2 and r = t\ + t 2 . If F(T) is an 
even function, then 



/ dh / dt 2 F(T) = - dT / drF(T) = 2 dT (t — T) F{T) . (33) 
Jo Jo 2 J-t J\t\ Jo 



With Eqs. (0), (0) and @, we may write Eq. (0]) as 

, fk(T 
(T - 2z) 



where 



and 



UT) = f y (T) = nHT + 2z y ( 35 ) 



8 a 2 (5 T 6 + 308 T 4 z 2 + 944 T 2 z 4 + 1728 z 6 ) (t - T) 
h[T) ~ 7r 4 (T + 2 z) 8 ' (36) 



The dispersion (: A-y 2 (í) :) in Eq. (|34] ) will be defined as a generalized principle value integral 
fT8| . Such integrals involving higher-order poles may be evaluated by successive integrations 
by parts, which remove the divergence at the point T = 2z and lead to the formula 



b 



p dx 



f(x) 1 



n-2 



-i b 



- 2 - i)! f^{x) (x - c )- n+1+i 

i=0 



(x-c) n (n-iy. 

+^TÏÏ^ fdxf^ix) (x-c)- 1 
(n — 1)1 Ja 



(37) 



We may now apply this formula to evaluate (: At> 2 (í) :). The result simplifies considerably 
if we assume that t ^> 2z, which is the limit of greatest physical interest. In this case, after 
some calculation one finds that 

(: ,):> = (: A^, ,) :) = ± (^ ^ + O (f ) (38) 

and 

/ a 2/ \ v 1 / 63 \ a 2 h 2 „ fz\ 

Note that even though there is no mean force in the x and y directions, the dispersion of 
the velocity in the direction parallel to the plate is still nonzero. 



S 



2. The cross term 



The other intriguing part of the quantum fluctuation of the Van der Waals force is the 
cross term. Its contribution to the velocity fluctuation, Eq. (|22|), is formally divergent. 
However, it can be made finite by an integration by parts procedure analogous to that used 
in the previous section. The key assumption which we need to introduce is one of adiabatic 
switching. This means that the effect of the plates is smoothly switched on in the past 
and then off in the future. Physically, this might be achieved by means of a plate whose 
reflectivity could be controlled. This switching will allow us to drop surface terms which 
would otherwise be divergent. An analogous switching was assumed in a treatment of the 



quantum fluctuations of radiation pressure ||15|| . There it was shown that the cross term plays 
a central role, and in fact gives the sole contribution when a làser beam in a coherent state is 
shined on a mirror. In this case, it was necessary to assume that the làser beam is switched 
on in the past and then off in the future in order to obtain finite velocity fluctuations for 
the mirror. 

In analogy to Eq. fl23D , the cross term of these force-force two point functions is defined 

by 



a 



(Fk(ti Z) F k (t , Z^ cross 



d k d k , <: Ei& : : E y E 3 ' :} £ 



(40) 



Due to Eq. (|Ï6|) , we need to know the vacuum two point function (EiEji) as well as the 
normal-ordered two point function (: E^Ej/ :) to compute the cross term. Use the equations 
(|) and (|7|) to compute these vacuum two-point functions Eq. (|Ï2|) , and find 



and 



(E x E x i)q 

(EyEy,) 

{E z E z i)q 



1 



2tt 2 
-1 



2vr 2 



-1 
2^2 



2 
2 

D* 



4(x-x') 2 4(í-í') 21 

D* + 
4(y - y" 2 



D 3 



+ 



D 3 
4(í-í') 2 
D* 



2 4(z - z'f 4(t-t') 2 
Jyi ^3 + ^3 



(41) 
(42) 

(43) 



Plug these vacuum two point functions, along with the normal-ordered two point functions 
Eqs. (|24|), (|25|) and ( pü|) into Eq. (p^6|) and compute the derivatives in Eq. (^). The 
force-force two point functions then becomes 



(F k (t,z)F k {t',z)) t 



a- 



Tí 



fi 



k.i 



4/, 



k.fi 



(44) 



Here f k ^ and f k ,6 are formed from the normal-ordered two point functions and their deriva- 
tives, and can be expressed as 



and as 



d k>k ,(: E x (x)E x ,{x') :) +d k u,{: E y (x)E y ,(x') :) +d k ^{: E z (x)E z ,(x') :) 



(45) 



x,6 



(: E x {x)E x ,{x') :) + 2(: E y {x)E y ,{x') :) + 2(: E z {x)E z ,{x') :) 



(46) 
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yfi 



and 



3.6 



2(: E x {x)E x ,{x') :) + (: E y {x)E y ,{x') :) + 2(: E z {x)E z ,{x') :) 



2(: ^.(x)^(x') :) + 2(: ^(x)^(x') :) + (: E z {x)E z ,{x') :) 



(47) 



(48) 



The singular parts, l/(íi — Í2) 4 and l/(íi — Í2) 6 , in Eq. (|44D are caused by the vacuum 
two point functions. Use Eq. ([44]) and change the variables (£1,^2) to dimensionless ones 
(si = ti/2z,S2 = t 2 /2z). The velocity fluctuation Eq. ([22] ) becomes 



(A«2(í)> 



cr 



cross o o 

7f z m z Jo Jo 



t/2z rt/2z 



fi 



k.\ 



4/ fc 



(2z) 2 (si - s 2 ) 4 {2 Z y{ Sl -s 2 y 



ds\ds2 



(49) 



Because of the adiabatic switching assumption discussed above, we can now integrate by 
parts and drop the surface terms, using the relations 



and 



W^W dsids2 = ~T2 



fk,6 i, 1 
:dSidS2 



{si - s 2 y 



240 



(d Sl ) 2 (d S2 ff k , 



(d sl nd S3 ) 3 f k , 6 



ln(si — S2) 2 ds\ds 2 



ln(si — S2) 2 dsids 2 



(50) 



(51) 



Plug these re-defined integrals into Eq. flUp and change variables to (u — s\— S2, v — S1 + S2). 
Use of the relation 



Jo 



ds\ds2 



\ f (0 ru+2a 

du 



a J —u 



dv + du 



a r2a—u 



Ju 



dv 



(52) 



leads to 

<Ai£(f)> 

where 



Q- 



rt/2z rt/2z 

/ / 9k(si,s 2 ) ln(si - s 2 ) 2 ds 1 ds 2 
Jo Jo 

/•í/22 rt/2z 

2t / gk(u 2 ) liíu 2 du — 2 / u gk(u 2 ) Inu 2 du 
Jo Jo 

3(5 + 275w 2 + 1325w 4 + 1041w 6 + 42w 8 ) 
16vr 2 (w 2 - lfz & 

3(23 + 663u 2 + 1573w 4 + 429u 6 ) 
: 8vr 2 (w 2 - 1) 9 ^ 8 ' 



cross n 

7T z m z 

a 2 



7T 2 m 2 



9x = 9y 



and 



9z 



(53) 



(54) 



(55) 



When t » 2z, the first term in Eq. ( [53"D goes to zero for all the components k = x,y,z, 
which leads to 



(Avl(z)) cross = (Av 2 (z)) a 
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13 h 2 a 2 
240 7r 4 m 2 z 8 



(56) 



and 

(a^)u ss - ■ ( 57 ) 

480 7r^m z z ò 

The magnitude of the velocity fluctuation of the z component is about 19 times of that of 
x and y components. In all cases, the contributions of the cross terms are larger than those 
of the fully normal-ordered terms, Eqs. (|38D and fl39|) 



(A^) cross = (Av 2 ) cross « 7.7 (: A^ : ) 



2/ 

(Av 2 z ) cross « -21 (: Ai; 2 :) . (58) 

The most surprising result is the negative z-component due to the cross term. The total 
velocity fluctuat ions are 



47 fi 2 a 2 

(A^ 2 ) = (Av 2 x ) = (: Av 2 x ■■) + (&v 2 x ) cross = (59) 



and 



<A# = <:A4:> + <Atf™ = -™£Hl. (60) 

The results are independent of time and the z component is still negative after the fully 
normal-ordered term is added to the cross term. The time-independent result shows a 
behavior similar to the case of Brownian motion in thermal equilibrium system. However 
we should also note that the velocity dispersion is not isotropic. The x and y components 
are much smaller than the magnitude of the z component 

(A«2> = (A«5>«0.06|(Atí>|. (61) 

The non-isotropic behavior is also reflected in the fact that the mean force is zero in the 
parallel direction, but non-zero in the perpendicular direction. That the (Av 2 ) approach 
constant vàlues, as opposed to growing in time, can be understood on the basis of energy 
conservat ion. 

Of particular interest is the fact that (Av 2 ) < 0. Recali that this quantity is a difference 
between a mean squared velocity with the plate and one without it, hence it is possible 
for this difference to be negative. (Similarly, the negative energy density in the Casimir 
effect arises from energy density being defined as a difference.) However, this negative value 
requires a physical interpretation. The most plausible explanation is that one cannot ignore 
the quantum nature of the test partides we have been discussing. The particle must have 
both a position uncertainty Az and a momentum uncertainty Ap z , obeying the uncertainty 
principle. Furthermore, because the particle is massive, there will be wavepacket spreading 
in which Az is an increasing function of time. Thus, even if the particle is initially in a 
minimum uncertainty wavepacket, at later time it will satisfy the uncertainty principle by a 
wide margin. Our interpretation of the negative (Av 2 ) is that the electromagnètic vacuum 
fluctuations cause a small reduction in the velocity spread of the wavepacket compare to 
what it would have been without the plate present. Imagine that we initially prepare the 
particle in a minimum uncertainty state, and then allow it to evolve for a time t » 2z). 
During this time Az will increase, but Av z will decrease slightly. 

In any case, the magnitude of the velocity changes due to electromagnètic vacuum fluc- 
tuations is always very small compared to the velocity spread due to quantum uncertainty. 
Let the latter be 
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Ap z Az . „ s 

Av q = — > — . 62 
m m 



Compare this to the spread due to vacuum fluctuations, 



Av f = max(v^A^) » . (63) 



Their ratio satisfies 



Av q V -2 / Vvr 2 z 3 

However, both factors on the right-hand-side of the above expression are small compared to 
one. The particle must be localized in a region small compared to the distance to the plate, 
so Az <C z. The size of the particle must also be small compared to z, and because the 
polarizability at is at most of the order of the volume of the particle, at <C z 3 . Thus 

^«1. (65) 
Av q v ; 



V. DISCUSSION 

Equations ( [59] ) and fl6"0|) teli us that the effect of the fluctuations of the retarded Van 
der Waals force is to generate a random motion around that described by the classical 
trajectory. Of course, if a particle is released in the vicinity of a conducting plate, it tends 
to fall toward the plate under the influence of the mean force, Eq. (0). However, we could 
apply a compensating classical force F c ; = — (F), so that the classical trajectory is that 
of a particle at fixed z. Nonetheless, it will still develop a mean squared velocity given 
by Eqs. fl38|) and fl39|). If we look at the x-direction (or y-direction), this is equivalent to 
thermal motion at a effective temperature of 

T e ff = i g , ( 66 ) 

V7687T 4 / k B mz* 

where k B is Boltzmann's constant. Equation flBlf ) can be written as 

— ^)(t) 8 (^ <-) 

where and a B are the mass and static polarizability of atòmic hydrogen, respectively. 
This effective temperature is essentially the temperature below which the system must be 
cooled so that the quantum fluctuation effects are not masked by ordinary thermal fluctu- 
ations. The effect in the z-direction is different from that in the transverse directions in 
that the mean squared velocity in that direction is reduced. Nonetheless, Eq. (|66D gives an 
estimate of the magnitude even in this case. The magnitude of the effect depends crucially 
upon how small z can be. For àtoms near a metal plate, both the assumptions of perfect 
conductivity and of using the static (as opposed to dynamic polarizability) break down for 
sufficiently small z, typically for z í$ 10 3 A. Thus the effect of the fluctuations will be very 
small in the range that both of these asumptions hold well. However, there is likely to be 
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some effect even at much smaller vàlues of z. A metal surface acts as a partial reflector of 
electromagnètic waves even up into the x-ray range, where Bragg scattering can produce 



reflectivities close to 100% at special angles Thus although Eq. flBSf) is strictly vàlid 
only for z > 10 3 A, it may produce crudely correct answers for z as small as as a few A. If so, 
the fluctuation effects could conceivably approach observable levels. This conjecture needs 
to be confirmed by more detailed treatements. 

The appearance of nonzero vàlues for (At^) and (At>^) requires some comment. By 
symmetry, a particle is equally likely to be deflected by the electromagnètic field in the +x 
or —x directions, and hence (v x ) = (v y ) = 0. However, the history of an individual particle 
does not have to respect the symmetry of the problem. Some partides acquire nonzero 
transverse components of velocity, leading to (At>^) 7^ 0. A similar situation arises in 
lightcone fluctuations due to quantum gravity effects in a compact space [ 2(J . Here Lorentz 



invariance holds on the average, but not for the history of an individual test particle. 

The time scale of the fluctuations due to the fully normal-ordered term are of the order of 
z, the light travel time between the particle and the plate, as may be seen from the fact that 
the correlation functions, Eqs. (|27| ) and (|28|), vanish for T ^> z. The time scale associated 
with the fluctuations arising from the cross term is of the same order. The short distance 
singularity of the cross term indicates that it contains fluctuations on arbitrarily short scales. 
However, these very ràpid fluctuations are averaged out by the time integrations. The final 
integral for (Avl(t)) cross , Eq. (^), again contains an integrand which vanishes rapidly for 
u = T/(2z) < 1. 

In summary, the Casimir-Polder result is a mean force, whereas the actual force is rapidly 
fluctuating. The typical magnitude of the fluctuations is of the same order as the mean force 
itself, but the time scale of the fluctuations is of the order of the light travel time between 
the atom and the plate. For most purposes, such as the Sukenik et al H experiment, 
the fluctuations average to zero and are not seen. In principle, it is possible to detect 
the fluctuations through the random motions which they will induce in test partides. For 
ordinary atòmic systems, this effect is very small. 

The effect discussed in this paper is also of interest in gravity theory. When quantum 
matter fields act as the souce of gravity, fluctuations of the stress tensor will lead to "passive" 
fluctuations of the spacetime geometry. These fluctuations are one of the physical phenomena 
to be expected in any quantum theory of gravity. Quantum fluctuations of the spacetime 
mètric imply Brownian motion of the test partides which probe the fluctuating mètric 
i3| , [i6|| . Thus the Brownian motion due to electromagnètic vacuum fluctuations treated 



here is a useful analogy for understanding the quantum nature of gravity. 
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